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ABSTRACT 

We present the first results of our HST WFPC2 F814W snapshot 
imaging survey, targeting virtually all sub-mJy decimetric radio-selected 
star-forming galaxies. The radio selection at ~ 1 GHz is free from ex- 
tinction effects and the radio luminosities are largely unaffected by AGN 
contamination, making these galaxies ideal tracers of the cosmic star for- 
mation history. A sub-sample of 4 targets is presented here, selected at 
1.4 GHz from the spectroscopically homogenous and complete samples of 
Benn et al. (1993) and Hopkins et al. (1999). The redshifts are confined 
to a narrow range around z ~ 0.2, to avoid differential evolution, with a 
radio luminosity close to L* where the galaxies dominate the comoving 
volume-averaged star formation rate. We find clearly disturbed morpholo- 
gies resembling those of ultraluminous infrared galaxies, indicating that 
galaxy interactions may be the dominant mechanism for triggering star 
formation at these epochs. The morphologies are also clearly different from 
coeval quasars and radiogalaxies, as found in star-forming galaxies se- 
lected at other wavelengths. This may prove challenging for models which 
propose direct causal links between AGN evolution and the cosmic star 
formation history at these epochs. The asymmetries are typically much 
larger than seen in the CFRS at similar redshifts, optical luminosities and 
Ha-derived star formation rates, indicating the possible existence of an 
obscuration-related morphological bias in such samples. 

Key words: cosmology: observations - galaxies: formation - infrared: 
galaxies - surveys - galaxies: evolution - galaxies: star-burst 



1 INTRODUCTION 

The study of the observational constraints on the cos- 
mic star formation history is currently among the 
most active fields in observational cosmology. The 
most widely used tracer of the comoving volume- 
averaged star formation rate (SFR) is the UV lumi- 
nosity density (e.g. Madau et al. 1996, Steidel et al. 
1996, 1999), found to peak at z ~ 1 — 2. However, 
little is known about the history of star formation in 
the Universe beyond its global average (e.g. Abraham 



et al. 1999). In this paper we present the first results 
from an on-going program with Hubble Space Tele- 
scope (HST) to study the morphology of starforming 
galaxies, using a sample of galaxies which dominate 
the comoving SFR and which were selected in a man- 
ner free from obscuration biases. 

An important caveat to the SFR constraints 
is that the UV luminosity of star-forming galaxies 
is dominated by the lowest-extinction regions. This 
leads to extreme sensitivity to obscuration corrections 
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Figure 1. HST F814W snapshots of the galaxies [CM84]144 (top left), [CM84]074 (top right), Phoenix-Deep-29 (bottom 
left) and Phoenix-Deep-96 (bottom right). Greyscale is linear, and contours are spaced in factors of 2. 
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[CM84]144 2 ' 3 ' 4 


08 55 47.60 


+17 02 28.4 
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18.7 


0.2253 


0.29 


18.3 


-21.2 


116.3 
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[CM84]074 2 > 3 ' 4 
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+17 03 46.9 


B = 


18.6 


0.2279 


0.29 


18.7 


-20.8 
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8.34 


0.112 


Phoenix-Deep-29 1 ' 5 ' 6 


01 10 43.88 


-45 51 22.7 


R = 


20.2 


0.210 


0.198 


19.7 


-19.6 


52.5 


1.1 


0.185 


Phoenix-Deep-96 1 ' 5 ' 6 


01 12 25.21 


-45 47 0.1 


R = 


18.4 


0.236 


0.225 


18.3 


-21.2 


75.9 


31.0 


0.317 



Table 2. HST z ~ 0.2 sub-mjy star-forming snapshot sample observed to date. All were observed for 2 X 400 seconds with 
the PC in filter F814W. References are (1) Georgakakis et al. 2000, (2) Condon k, Mitchell 1984, (3) Rowan-Robinson et 
al. 1993, (4) Benn et al. 1993, (5) Hopkins et al. 1998, (6) Hopkins et al. 1999. Star formation rates (1.4 GHz and Ha) are 
in Mq yr — 1 , and A are the asymmetries. For comparison, the z = 0.2 (z = 0.3) radio L t assuming (1 + z) 3 luminosity 
evolution corresponds to 63.4 Mq yr _1 (80.6 Mq yr _1 ). The Tresse &; Maddox (1997) Ha L* at z ~ 0.2 corresponds to 
22.4 M yr" 1 . 
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Table 1. Flux of the break in the 1.4 GHz luminosity 
function (L*) under various luminosity evolution assump- 
tions. Uses the luminosity function from Serjeant et al. 
1998, with an n = 1, A = world model. Very similar 
results are derived from the Condon 1992 luminosity func- 
tion. For a Hubble constant of Ho = 50 km s — 1 Mpc — 1 
the characteristic luminosity lies at ~ 2.8 X 10 22 W Hz -1 . 
The fluxes themselves are Ho-indcpcndcnt. 

(e.g. Pettini et al. 1998, Meurer et al. 1997), large 
enough to eliminate the evidence for a redshift cut-off 
in the SFR. Indeed, SCUBA observations of the Hub- 
ble Deep Field (HDF) (Hughes et al. 1998) detected 
5 ultraluminous galaxies at z > 2 (the redshift con- 
straints come mainly from the radio-FIR correlation 



and not from the uncertain HDF IDs), which together 
implied an obscured star formation rate at least as 
large as the integrated de-reddened UV-derived rate. 
Moreover, the Ha-derived star formation rate from 
the Canada-France Redshift Survey (CFRS, Tresse & 
Maddox 1997) was a factor ~ x 3 larger than the UV 
estimate in the same sample, but comparable to those 
derived from the mid-IR for the CFRS (Flores et al. 
1999) or the HDF (Rowan-Robinson et al. 1997). 

However, a serious problem is that optical-UV 
light is strongly skewed to low extinction regions, 
while the reverse is true for mid and far-IR selected 
samples, leading to problematic obscuration-related 
effects in either case. An unbiased technique for se- 
lecting star-forming galaxies is to sample at decimet- 
ric radio wavelengths, where both obscured and un- 
obscured star forming galaxies contribute (Condon 
1992). Several groups (including ourselves) have be- 
gun to exploit this technique to trace the comic star 
formation history, and to study the obscuration ef- 
fects in star forming galaxies (e.g. Cram et al. 1998, 
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Figure 2. Asymmetry parameters for the radio-selected 
sub-mJy star-forming galaxies (filled symbols) compared 
to the CFRS/LDSS galaxies with 0.17 < z < 0.3 (open 
symbols). 



Serjeant et al. 1998, Oliver et al. 1998, Haarsma & 
Partridge 1998, Cram 1998, etc). 

We have embarked on a programme to image 
almost all star-forming galaxies selected at decimet- 
ric radio wavelengths, using the WFPC2 in snapshot 
mode on the HST. In this paper we carry out a mor- 
phological study of a small subset of these galaxies 
at z ~ 0.2, from our first cycle 8 observations in this 
programme. (Although these are at low redshifts by 
many standards, they are not local, with the volume- 
averaged SFR being a factor of ~ x2 higher than in 
the local Universe for (1 + z) 3 luminosity evolution 
models (Table [j]).) The galaxies in this sample are se- 
lected so that they dominate the (radio-derived) star 
formation history of the Universe at this epoch, i.e. 
close to the L* characteristic luminosity in the radio 
luminosity function. 

The sample selection is discussed in Section 2, 
which also presents the HST observations and data 
reduction. These data are then analysed in Section 3, 
followed by a discussion of the results in Section 4. 



2 SAMPLE SELECTION AND 
OBSERVATIONS 

The following criteria are considered in selecting the 
sample for HST observations: 

(i) The objects are sub-mJy radio sources, selected 
at 1.4 GHz (Benn et al (1993); Georgakakis et al 
(2000)). This avoids obscuration-related selection bi- 
ases. Decimetric radio fluxes are much less affected by 
AGN contribution than shorter wavelength observa- 
tions (compare e.g. the ongoing contoversy over the 
AGN fraction in higher frequency sample with Ham- 
mer et al. (1995) finding ~ 50% AGN fraction at 5-8 
GHz while Windhorst et al. (1995) find only ~ 5%). 

(ii) Spectroscopic data are used to select emission- 
line star-forming galaxies. 



(iii) The galaxies must not lie close to bright stars 
to avoid HST roll angle constraints. 

(iv) Galaxies selected to lie close to the 1.4GHz L, , 
which dominate the volume averag ed SFR (Table [TJ> . 

(v) Galaxies in a narrow redshift range around z ~ 
0.2 to avoid differential evolution. 

Our primary sample satisfies (i), (ii) and (iii), and 
is the subject of our 150 snapshots in cycles 8 and 9. 
Ten galaxies also satisfy conditions (iv) and (v); these 
are the cycle 8 targets. The HST observations for four 
of these have now been completed. We selected the 
F814W (wide I band) filter. At this wavelength, the 
light is mostly dominated by old stellar population, so 
gives a more accurate measure of distortion in the un- 
derlying gravitational potential than e.g. UV observa- 
tions. Details of the current sub-sample are presented 
in Table §. 

The data were reduced using the Interactive Data 
Language (IDL), starting from the automatic pipeline 
products. Cosmic rays were identified as > 3<j differ- 
ences between frames. For sky subtraction, we esti- 
mated the modal value of the underlying counts distri- 
bution using iterative fits to the readout histograms. 
The final HST images for the four sub-mJy radio 
sources are presented in Figure |l| 



3 RESULTS 

[CM84] 144 is extremely disturbed, with a tidal tail 
extending ~ 1.5" North-East of the nucleus. Sev- 
eral secondary nuclei are also apparent, reminiscent 
of HST WFPC2 F814W imaging of ultra-luminous 
infrared galaxies (e.g. Borne et al. 2000) such as Arp 
220 (Borne & Lucas 1997). 

[CM84] 074 is only mildly asymmetric. There are 
hints of face-on spiral structure, and an excess flux 
~ 0.5" West of the nucleus. There are two companion 
galaxies a few arcseconds to the West, also visible in 
the digitised sky survey. (We reject these as identifi- 
cations of the radio source based on the ID magnitude 
quoted in Benn et al. 1993.) 

Galaxy Phoenix-Deep-29 is again not as clearly 
disturbed. However, the inner isophotes are clearly 
offset in position and (tentatively) in orientation from 
the outer isophotes. This galaxy also shows signs of 
interaction with a companion 5" to the East. There 
is a clear tidal tail from the companion extending a 
few arcseconds North. 

Galaxy Phoenix-Deep-96 shows clear signs of 
morphological disturbance, with the structure dom- 
inated by a bright central bar-like feature. There are 
also hints of multiple nuclei inside the bar, and there 
is a secondary peak ~ 0.5" South- West of the bar. 
There are however no other clearly associated com- 
panion galaxies visible in the HST image. 

To quantify the morphological disturbance, we 
follow the procedure adopted by Brinchmann et al. 
(1998) and others, by rotating each image through 
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tt and subtracting it from the original image. Nor- 
malising this to the galaxy flux yields the fractional 
asymmetric flux: 



.4 



£|G, 



EG, 



k 



(1) 



where G (G') is the rotated (unrotated) image, k is 
a correction for systematics such as sky gradients, 
and the sum is performed over pixels ij. The values 
of A are fairly insensitive to the choice of aperture 
used for the self-subtraction, provided all pixels at 
> 1.5cr are included in the subimage. The magnitude 
of k can be estimated by applying the technique to 
blank sky regions of the same size. The errors in A 
are in practice dominated by the uncertainty in k: we 
find typically AA ~ 0.05. A further advantage of this 
statistic is that the /-selected Canada-France Red- 
shift Survey (CFRS) and bj-selected LDSS surveys 
have extensive published HST morphologies quanti- 
fied using the same statistic, yielding valuable control 
samples. The asymmetry measures for the sub-mjy 
star-forming galaxies, quantified as discussed above, 
are listed in Table El 



4 DISCUSSION AND CONCLUSIONS 

In Figure ^| we plot the asymmetry statistics of our 
galaxies, and compare them to a control sample of 
coeval galaxies from the CFRS sample (Brinchmann 
et al. 1998). The control extends to fainter 7fsi4w 
fluxes than our star-forming galaxy sample, though 
there is some overlap. Remarkably, our galaxies are 
far more asymmetric than the control sample of an 
(ostensibly) identical population of optically-selected 
star- forming galaxies of the same Jfsi4W- 

The optical control must be well-matched with 
the radio sample in as many optical properties as 
possible. The radio sample is more asymmetric at 
a fixed optical luminosity, but is this also true at a 
fixed optically-derived SFR? I.e., could the morpho- 
logical differences between the radio sample and the 
optically-selected control be instead because they are 
not well-matched in optically-derived star formation 
rates? To address this we must compare the Ha lumi- 
nosities of our radio sample with those of the control. 
The data for the radio samples are presented in Benn 
et al (1993) and Georgakakis et al (2000). The indi- 
vidual Ha luminosities of the CFRS galaxies are not 
published, but Tresse & Maddox (1998) report a tight 
correlation between Ha luminosities and Ms (AB) ab- 
solute magnitudes. In Figure ^ we use this relation 
to estimate the Ha-derived star formation rates, and 
compare the optical control sample with our sub-mjy 
star-forming galaxies. We adopt the following con- 
version between Ha luminosities and star-formation 
rates (SFRs): 



which is derived assuming Salpeter IMF between 0.1 
M a and 125 M (see Serjeant et al. 1998, Oliver et 
al. 1998). We de-redden our Ha fluxes by Av = 1 for 
consistency with optically-selected samples. (Tresse & 
Maddox 1998 note that the Ha luminosity function is 
essentially unchanged if assuming this Av through- 
out their sample, instead of using individual Balmer 
decrements.) There is a more substantial overlap of 
the Ha-derived SFR with the control sample, and 
it is again clear that at a fixed Ha SFR the radio- 
selected sample is significantly more disturbed than 
the optically-selected sample of star forming galaxies. 

In summary, the radio-selected sources are typ- 
ically more morphologically disturbed than the 
optically-selected control, well-matched in redshift, 
optical luminosity and/or Ha SFR. Why this might 
be the case is hinted in Table ^, where we compare the 
star formation rates of our galaxies estimated from 
Ha and the 1.4GHz radio luminosity. To convert deci- 
metric radio luminosities to SFRs we use 



SFR(1.4GHz) = Li. 4 /(7.63 x 10 20 W Hz" 1 ) 



(3) 



appropriate for the same IMF as above (see Oliver 
et al. 1998, Serjeant et al. 1998). Our Ha SFRs are 
comparable in most cases to L* in the Ha luminosity 
function (Tresse & Maddox 1997). However the radio- 
derived star formation rates are much larger (though 
still comparable with the radio L», Table hi). This in- 
dicates that a large fraction of the star formation in 
these objects occurs in regions with Av 3> 1, since 
the observed Ha luminosities will be dominated by 
regions with Ay ~ 1. Indeed if the dust is well-mixed 
with the Ha emitting gas, the different optical depths 
for Ha and H/3 ensure that Av = 1.1 would be de- 
rived for a simple screen Balmer decrement model 
ss of the true extinction to the rear of the 



SFR(Ha) = L(Ha)/(1.41 x 10 34 W) 



(2) 



cloud. We therefore suggest that optically-selected 
samples may under-represent disturbed galaxies with 
large amounts of obscured star formation, in any SFR- 
weighted quantity. 

The galaxy morphologies of the sub-mjy sources 
are also markedly different from coeval AGN. McLure 
el al. (1999) and McLeod & Reike (1995) have found 
giant ellipticals hosting both radio-loud and radio- 
quiet quasars at these redshifts, and radiogalaxies are 
also early type galaxies at these epochs. The EMSS- 
selected sample of Schade et al. (1999) at z < 0.15 
also finds no evidence for strong interaction/merger 
activity for any AGN in the sample. This may pose 
problems for models which propose causal links be- 
tween the cosmic evolution of AGN and the cosmic 
star formation history. For example, a model in which 
both are driven by similar merger events must also 
provide a mechanism for delaying the onset of central 
engine fuelling, to account for the (relatively) more re- 
laxed AGN host galaxy potentials. This is related to 
the more fundamental and unsolved problem of AGN 
evolution: how to drive the gas fuel down ~ 5 orders 
of magnitude in radius. 
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Figure 3. Asymmetry parameter plotted against the es- 
timate of the Ha star formation rate discussed in the text. 
Symbols as in Figure ^| 



We can also compare our sample with the ultralu- 
mirious galaxy samples of e.g. Borne et al. (2000), who 
also used HST WFPC2 F814W in snapshot mode. 
These authors found highly disturbed systems often 
in moderately rich environments, which they used to 
argue for an evolutionary sequence involving com- 
pact groups and ultraluminous galaxies. Although 
asymmetry statistics have not been calculated for the 
ultraluminous galaxies, our targets appear qualita- 
tively less disturbed and the environments possibly 
less rich. (Ideally both ultraluminous galaxies and 
radio-selected star-forming galaxies should also have 
their environments quantified with for example the 
B gg statistic; e.g. Yee & Green 1984, Hill & Lilly 1991, 
Wold et al. 2000) This may suggest a link between 
the level of star formation activity after interactions 
and the richness of the environment, perhaps via the 
number or the nature of the interaction events. 

Although our current sample is small, the asym- 
metric morphologies in our ~ L, targets suggest that 
galaxy interactions play a major role in the evolution 
of the star formation and metal production rates in 
the low-z Universe. Further papers in this series will 
present results for larger samples of sub-mJy star- 
forming galaxies. However, it is worth keeping in mind 
that the formation and evolution of galaxies cannot 
be characterised by a single parameter, such as the 
volume- averaged star formation rate. In expanding on 
this simple first-order description, the obscuration- 
independent selection of galaxies will be essential, 
as will co-ordinated multi-wavelength follow-up (e.g. 
ELAIS, Oliver et al. 2000). 
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